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IN MEMORY OF OKKIE DE JAGER 

In a Pulsar Wind Nebula (PWN), the lifetime of inverse Compton emitting electrons exceeds the 
lifetime of its progenitor pulsar, but it exceeds also the age of the electrons that emit via synchrotron 
radiation; i.e. during the evolution of the PWN, it can remain bright in IC, whereas its GeV- TeV 
gamma-ray (for 10^-10^ years) flux remains high for timescales much larger than the Pulsar lifetime 
and the PWN visible in X-rays. The shell-type remnant of the supernova explosion in which the 
pulsar was formed also has a much shorter lifetime. In this scenario, the magnetic field in the cavity 
induced by the wind of the progenitor star plays a crucial role. This is in line with the discovery 
of several unidentified sources in the TeV gamma-ray band without X-ray counterparts. Moreover, 
the consequences are important also in order to reinterprete the detection of starburst galaxies in 
the TeV gamma-ray band considering a leptonic origin of the gamma-ray signal. 



I. PULSAR WIND NEBULAE 

Pulsar wind nebulae (PWN) represent a unique lab- 
oratory for the highly efficient acceleration of particles 
to ultra-relativistic energies. Early ideas about parti- 
cle acceleration by electromagnetic waves emitted by 
pulsars due to [i| led to the discovery by [2|] that a 
significant fraction of the spin-down power of a pulsar 
such as Crab is dissipated through a relativistic wind 
with in situ acceleration of particles at its termination 
shock, since the cooling time scales of the synchrotron- 
emitting electrons are very short in X-rays. Further 
developments led to the model of [s"] involving a low 
level of magnetization at the shock, raising the ques- 
tion of the transition from high to low sigma since the 
launch of the wind consisting of pairs and Poynting 
flux involves a high magnetization 
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where p denotes the density and 7 the Lorentz fac- 
tor of the wind electrons and positrons. Recently, 
MHD simulations of such winds were performed by 
[H [5| . Central of our current understanding of PWN 
is that they are equatorial and highly relativistic, 
showing features such as backflows and jets. This 
might provide clues to their observational identifica- 
tion, since energy release presumably due to recon- 
nection can lead to flares from the wind, considering 
relativistic transport effects and beaming, as recently 
observed in Crab. Moreover, the MHD models provide 
clues for predicting the temporal evolution of PWN, 
and to bring this into accordance with the observa- 
tions of a large sample of putative PWN of different 
ages. 



II. ANCIENT PULSAR WIND NEBULAE 

Although they are often detected as non-thermal 
X-ray sources, an evolved PWN can indeed lead to 
a fairly bright 7-ray source without any lower en- 
ergy counterpart @ Q. The key issue is that the 
low energy synchrotron emission, where te is the syn- 
chrotron emitting lifetime of TeV 7-ray emitting lep- 
tons with energy E, depends on the internal PWN 
magnetic fleld [3j which may vary as a function of 
time, following te (x t^" if B{t) cx where a 

is the power-law index of the decay of the average 
nebular field strength, whereas the VHE emission de- 
pends on the CMB radiation field, which is constant 
on timescales relevant for PWN evolution. MHD sim- 
ulations of composite SNRs [1] find that a = 1.3 until 
the passage of the reverse shock; after it is expected 
that field decay would continue for much longer time 
since expansion continues even after the passage of 
the reverse shock. From hydrodynamic simulations, 
an expression for the time of the reverse shock pas- 
sage (i.e. the return time of the reverse shock to the 
origin) has been given Q as 
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where the first term is the density of the ISM, the 
second the ejecta mass during the SNR explosion and 
the third the SNR blast wave energy. The stellar wind 
of a high-mass star can blow a cavity around the pro- 
genitor star [1^ with relatively low ISM density, so 
that Tfj >> 10 kyr. In such a case it is expected that 
B{t) can decay as until the field is low enough 
for the X-ray flux to drop below the typical sensi- 
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FIG. 1: Representation of the Spectral Energy Distribu- 
tion of G21.5-0.9 during the first 25 kyrs from [12]; this 
example describes efficiently the temporal behavior of syn- 
chrotron and IC components 



microgauss in the inner Galactic disc, and presum- 
ably drops to tenths of microgauss at higher latitudes 
[U [H]. Although neglected in previous studies, this 
magnetic field of the ISM into which the Pulsar Wind 
Nebula is expanding must have some effect onto the 
morphology and time evolution of the PWN. As the 
magnetic field in the PWN drops rather rapidly over 
time, as , it becomes comparable to the sur- 

rounding ISM field at timescales roughly compara- 
ble to the transition timescale into the post-Sedov 
phase. However, as noted above, different PWN will 
be located in environments where the ambient field 
strength is different by factors of order unity, depend- 
ing rather strictly on the location of the PWN in the 
Galaxy. The respective timescales for the field to be- 
come comparable to the ISM field will therefore differ 
by a factor of ^ 3 for a PWN at the inner disc com- 
pared to one in the Galactic halo. Detailed studies 
of the effects on PWN morphology and evolution are 
highly demanding and yet challenging, owing to the 
limited resolution of high energy instruments and the 
relative scarcity of identified PWN. 



IV. TWO OTHER IMPORTANT 
COROLLARIES 

A. TeV unidentified sources 



tivity levels. As a result, in a scenario where the 
magnetic field decays as a function of time, the syn- 
chrotron emission will also fade as the PWN evolves. 
The reduced synchrotron losses for high-energy elec- 
trons for such a scenario will then lead to increased 
lifetimes for these leptonic particles. For timescales 
shorter than the inverse- Compton lifetime of the elec- 
trons {tic oc 1.2 X l{)^{Ee/neY)-^ years), this will 
result in an accumulation of VHE electrons which will 
also lead to an increased gamma-ray production due 
to up-scattering of CMB photons. Such accumulation 
of very-high energy electrons in a PWN has indeed 
been seen in many TeV PWNe (such as HESS J1825- 
137 [HI). To summarize, during their evolution PWN 
may appear as gamma-ray sources with only very faint 
low-energy counterparts and this may represent a vi- 
able model for many unidentified TeV sources. This 
effect can be clearly seen in Fig. [H shown by Okkie 
de Jager in SciNeGHE 2008 conference 



III. LINK BETWEEN PWN TIMESCALES 
AND GALACTIC LATITUDES 

There are numerous observational hints pointing 
to the fact that the interstellar magnetic field in the 
Milky Way has field strength of the order of several 



Ancient PWNe are a natural explanation for TeV 
unidentified sources Q . In fact ~ 50% of TeV Galac- 
tic sources are still formally unidentified |1 5|; almost 
all (among them, only HESS J0632-|-057ll6ll, that 
is considered a 7-ray binary candidate, and HESS 
J1943+213 [l3|, considered a HBL candidate but re- 
cently suggested to be more likely a PWN [1^, are 
point-like) are extended objects with angular sizes 
ranging from approximately 3 to 18 arc minutes (but 
of few more extended exceptions, such as HESS J1841- 
055 d^), lying close to the Galactic plane (suggest- 
ing a location within the Galaxy). In each case, the 
spectrum of the sources in the TeV energy range can 
be characterized as a power-law with a differential 
spectral index in the range 2.1 to 2.5. The gen- 
eral characteristics of these sources (spectra, size, and 
position) are similar to previously identified galactic 
VHE sources (e.g. PWNe), however these sources 
have so far no clear counterpart in lower-energy wave- 
bands. This scenario is particularly suitable for sev- 
eral sources that, even after deep MWL observa- 
tions, lack any plausible X-ray counterparts (such as 
HESS J1507-622 [?], HESS J1427-608^ESS J1708- 
410 ^] and HESS J1616-508 [l^l (Hi). On the 
other hand, several unidentified sources have already 
been identified as PWNe after their discovery (such as 
HESS J1857+026 [Hi or HESS J1303-631 [22]) and in 
sources that have several plausible counterparts, the 
PWNe contribution can hardly be avoided (such as 
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FIG. 2: Preliminary representation of the Spectral En- 
ergy Distribution of HESS J1507-622 assuming initial con- 
ditions similar to the G21.5-0.9/PSR J1833-1034 system; 
the radio information is taken from [SJ, the Fermi-LAT 
points from [2^, the Chandra and H.E.S.S. points from 



hot spot B in HESS J1745-303 23] or HESS J1841- 
055). 

Fig. E] shows the example of HESS J1507-622: this 
VHE source does not have any SNR as a possible coun- 
terpart and does not show any coincident or close- 
by pulsar [3]; as said, it is expected that the pulsar 
already spun-down and that the X-ray nebula faded 
away below the sensitivity of current X-ray instru- 
ments, however it requires an arbitrary choice of the 
initial conditions in order to model this source: in the 
shown example we assume to have conditions similar 
to the G21.5-0.9/PSR J1833-1034 system. 



B. Starburst galaxies 

Another important consequence of these long-living 
gamma-ray sources regards starburst galaxies: in fact 
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FIG. 3: This figure taken from [2^ shows the predicted 
TeV gamma-ray luminosity from PWN (sohd line) for dif- 
ferent core-collapse supernova rates. Note that the points 
shown for LIVIC and the Milky Way are extrapolations of 
measurements obtained at Fermi energies. 



it has been recently shown [2^ that PWNe are impor- 
tant and not negligible in explaining the TeV emission 
detected from NGC 253 [13 and from M82 [28]. In 
fact, even if M82 and NGC253 show a much harder 
GeV gamma-ray spectrum than the Milky Way (in 
line with the high gas density in the star forming re- 
gions (SFRs) and a cosmic-ray origin of the gamma 
rays), diffusive-convective escape of the cosmic rays 
from the SFRs should lead to a steepening of the 
cosmic-ray induced gamma-ray spectrum above ~ 10 
GeV, in which case the cosmic rays would fall short in 
explaining the TeV luminosities; on the other hand 
PWNe associated with core-collapse supernovae in 
SFRs can readily explain the observed high TeV lu- 
minosities. The proof of this could arrive from deeper 
gamma-ray observations on other galaxies, as shown 
in Fig. H 

Moreover also neutrino observations could reveal if 
the PWN play a role in accelerating cosmic ray pro- 
tons and ions; the highly magnetized pulsar equato- 
rial wind may also give rise to non-negligible accel- 
eration of cosmic rays by this source class. Protons 
and ions, swept into the plerion by the reverse shock, 
could be accelerated efficiently at the ultrarelativistic 
pulsar wind shocks ^29,] ^] . This scenario could pos- 
sibly be tested by searching for high-energy neutrinos 
from galactic PWNs like the Crab. Until now, such 
measurements have turned up only null-results. How- 
ever, a stacking analysis of candidate sources could 
significantly enhance the neutrino telescope sensitiv- 
ity, thus resulting in a much more decisive test of this 
scenario. Moreover recently strong upper limits have 
been given by means of Ice Cube (e.g. jlH), limits 
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close to the theoretical predictions. 
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